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Exact and explicit probability densities for one-sided Levy stable distributions 
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We study functions gd^) which are one-sided, heavy-tailed Levy stable probability distributions 
of index a, < a < 1, of fundamental importance in random systems, for anomalous diffusion 
and fractional kinetics. We furnish exact and explicit expressions for gdx), < x < oo, for all 
a — l/k < 1, with k and I positive integers. We reproduce all the known results given by fc < 4 
and present many new exact solutions for fc > 4, all expressed in terms of known functions. This 
will allow a 'fine-tuning' of a in order to adapt ga{x) to a given experimental situation. 
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Theoretical description of many collective physical sys- 
tems which include special sort of disorder or randomness 
often requires a radical departure from classical diffu- 
sive behaviour. On the probabilistic level this signifies 
the appearance of distributions with non-conventional 
characteristics, like diverging mean and variance along 
with all integer moments different from zeroth one. 
In this context the prominent role plays the discovery 
of particular distributions with these properties, called 
now Levy stable laws [l[, whose generic example is 
51/2(2;) = (2i/7rx^/2)"i exp(-l/4a;), x > 0; the word 
"stable" means here that the product of characteristic 
functions (cf) of two such laws is a cf of another law 
of the same type [1|. The general distribution of that 
type ga {x) can be shown to possess the cf or the Laplace 
transform of the form & 



e ^^ ga{x) dx 



p>0,0<a<l, (1) 



which is the well-known Kohlrausch- Williams- Watts 
function [5| or stretched exponential. Several indepen- 
dent proofs can be given that ga{x) obeying Eq. ([1]) is 
positive UMM- 

The functions ga{x) are ubiquitous in many fields of 
condensed and soft matter physics [3-Q , geophysics llOll . 
meteorology [11|, economics [i2l, fractional kinetics [Ij, 
[l4| . etc. For instance, the value a = 1/4 is thought to 
describe mechanical and dielectric properties of glassy 
polymers [la] ■ It is also confirmed that the same value of 
a is relevant for a statistical description of subrecoil laser 
cooling [ly, ll7| . In general, numerous phenomena falling 
in the class of subdifFusion [1^ call for ga{x), a < 1, 
in their theoretical description. On theoretical side, the 
Levy stable distributions are essential tools in the study 
of random maps and resulting combinatorial structures 
[19| . The actual use of Levy stable type distributions has 



been hampered for subjective and objective reasons [20, 
[2]| . The subjective ones include a certain reticence to 
use distributions with both mean and variance diverging. 
The main objective reason is a lack of knowledge of 170(0;) 
for most values of a. The existing interpolation formulae 
[iSj appear to be cumbersome to use. 

It seems that obtaining explicit ga{x) for arbitrary 
< a < 1 constitutes a true challenge: only for a 
limited number of values of a, i.e. a = 1/2 (s. 51/2(2;) 
above), 1/4 [H, 1/3 [H, 2/3 [H, and 3/4 ^ the ex- 
plicit forms of Qa (x) are known. The formal solution for 
arbitrary a [j, |9| is only of a limited use as it requires 
series or asymptotic expansions, which may become prob- 
lematic, especially for small a. 

The objective of this work is to present an universal 
formula for ga{x), a = l/k, with k > / positive integers, 
which is exact and explicit. It reproduces all the known 
cases enumerated above, and yields an infinity of new 
solutions for fc > 4, of which we quote, for the first time, 
several instances. 

Eq. ([1} for a = l/k can be inverted giving 



91/ k[ 



(27r)('=^ 



0/2 ^ ^' 



'Lk 



I' 



k^ : 



A((,0) 
A(fc,0) 



(2) 



vahd for all a; > 0, where G^"^" (^I(m) i^ the Meijer G 

function [H, [3^ and A(fc, a) = f , ^, . . . , ^i±|^ is a 
special list of fc elements. Eq. ^ is listed without proof 
as a special case for z/ = and a = 1 of formula 2.2.1.19, 
in vol. 5 of [30|- The detailed demonstration of Eq. ([2]) 
with a combined use of Laplace and Mellin transforms 
will be given elsewhere. It turns out that the r. h. s. of 
Eq. ((2|) is a finite sum of fc — 1 generalized hypergeometric 

functions of type pFq i Z^"} | z j [30] : 
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9i/k{x) 






b j{k,l ) 



1+ 



iFk 



l,A{l,l+jl/k) 
A{k,j+1) 



(-1) 



k-l 



I' 



kf": 



(3) 



where bj{k,l) are numerical coefficients given by 



b,{k,l) 



P'/'= vii 



k^ (27r)(fc-0/2 



UlZlTC-^) 




nt-^ir(^); 



ntlr(i-i) 



(4) 

where r(?;) is Euler's gamma function. Eq. ^ is the ex- 
act implementation of the program outlined in the fun- 
damental work of Scher and Montroll [231, i^ which it 
was conjectured thai Qj/k{x) can be expressed in terms 
of p-Fq's. However in [23| actually only one new instance 
53/4(2;) was written down. Our formula Eq. ([3]), after ap- 
propriate reductions in pFg's, see below, gives all exactly 
known cases mentioned above [20l - |23t . with 51/4(0;) |22| 



51/4(2;) 



&i(4,l) 



,5/4 



0-^2 1 1/2,3/4 



-1 



&2(4,1) 



^3/2 



(5) 



X 0^2 



2 I 3/4,5/4 



-1 



&3(4,1) 



.F; 



7/4 0-^2 I 5/4,3/2 



-1 



and offers an unlimited number of new solutions for 
9i/k{x), fc > 4, e.g.: 



5p/5(a;) 



E 



„l+jp/5 P+ 



iJ^E 



A,A(p,l+jp/5) 
I A(5,j + 1) 



P^ 



^^xP J 



,(6) 



p — 1, . . . , 4, see Table U for coefficients in Eqs. ([5]) and 
®, etc. 

The symbol A(A;, a) in Eq. ([3]) permits one to en- 
code all the possible cases of k and I in a single for- 
mula. However we draw attention to the fact that can- 
cellations will appear there due to the obvious identity 



p+r-f*g+ 



/(op 



),(«r.) 



p-Tq 



/(Oj 



where (a^) is an 

arbitrary sequence of r parameters not equal to zero or 
to negative integers. Thus Eq. ([5]) is a sum of three 
qF2 functions, and likewise 53/4(2;), which is not spec- 
ified here, will be a sum of three 2^2 functions, neatly 
eon&ming Eq. (CIO) of [23, etc. In this manner for any 
Z/fc, a closed form oi gi/k(x) can be obtained from Eq. ([3|). 
However, only for fc < 3 it can be written down in terms 
of standard special functions 0, [H [HI, [H, [13 . 

A heuristic indication of how Eq. ([3]) comes about can 
be obtained from the series representation for ga{x) de- 
rived by Humbert |25l|. discussed for example by Hughes 
1291 and used in [13 ■ The series 



la{x) 



i V 



Uj'-^ 



(■_iy+i 

^^^r(H-aj)sin(7r«j-) (7) 



is a convergent expansion valid for all < a < 1 and 
X > 0. For a = l/k the decomposition of the summation 
index j modulo fc yields an equivalent representation of 
gi/k{x) as a sum of fc — 1 inffiiite series. The structure 
of the coefficients in each of these series involve gamma 
functions ratios of type T{li + 9)/r{ki + (j)), where i is the 
new summation index and 6 and (f) are simple functions 



of fc and /. The application of Gauss-Legendre multipli- 
cation formula to both of these gamma functions allows 
the identification of pFq^s in Eq. ([3]) and the extraction 
of coefficients bjik.l) of Eq. (ji]). 

The advantage of our solution, Eqs. ([3]) and (jj]) over 
Eq. ([7]) is clearly seen in practice, in conjunction with 
the use of computer algebra systems J27| . Since in re- 
cent versions of these systems the hypergeometric func- 
tions pFq as well as the Meijer G function are fully im- 
plemented, their use permits high-precision calculations. 
For reader's convenience we give in [2^1 the Maple® syn- 
tax for gi/k{x), see Eq. © above. Our experience indi- 
cates that for small a our results for small x are more 
practical to use than the a; — )■ asymptotics given in 
(3. The reason is that there the region of applicability 
of Mikusinski's asymptotic expansion |3| shrinks to ex- 
ceedingly small values of x. For example for a = 1/20, 
5i/2o(0) = 0, but a huge peak in 51/20(2;) appears already 
at a; ^ 10~^4, jn contrast, our formulae work fine for any 
X in this region. In the opposite limit for a < 1 the Hum- 
bert expansion Eq. ([7]) is slowly convergent for x < a but 
approximation [3| works well as then ga{x) is very close 
to zero in a considerable region near x = (e.g. already 
for a = 5/6 the function 55/5(2;) is practically equal to 
zero up to a; ~ 0.35). Such a practically flat region for 
small x can also be seen for a ~ 4/5, compare the curve 
HI on Fig. [31 

In Fig. [T] we compare three distributions for l/k = 1/2, 
1/3 and 1/4. The salient feature for l/k = 1/4 is the 
appearance of a sharp maximum for very small x so that 
these three curves can be barely shown on the same scale. 
Analogously, for l/k ~ 1/5 the maximum of 51/5(2;) ap- 
pears at a;o(l/5) « 0.0002 and the value 51/5(2:0) ~ 25. 
For X < a;o(l/5) the values of 51/5(2;) are very close to 
zero. As already mentioned above, for smaller values 
of l/k this type of behaviour is even more pronounced 
and it explains a posteriori the difficulties encountered 
in devising approximations valid for small l/k and small 
X [la, [23l- In Fig- ID we present the comparison of several 
distributions for values l/k « 1/2. Here the 'sharpen- 
ing' of the distributions, as l/k goes from 1/2 to smaller 
values, is very clearly visible but is less dramatic than in 
Fig.[T] We present in Fig.[3]the new distributions 5p/5(2;) 
given by Eq. ([6]) for p = 2, 3 and 4. 

All these probability distributions share the following 
features: (a) 5a (a;) — >■ 0, for a; — >■ 0, where they present 

-2+a -a 

an essential singularity ~ a;^(i-°> exp[— A(a)xi-°1, 
A(a) > [3; (b) g^ix) -^ B{a)x-^^+°'\ for 
X -^ 00, B(a) > 0, indicating heavy-tailed asymp- 
totics for large x; (c) all their fractional moments 
Ma(/i) = J^ x''gaix)dx = r(-^/a)/[ar(-/x)l, for 
real /i, —00 < fj, < a, including Ma{0) = 1, are finite, 
and are infinite otherwise; (d) ga{x) are unimodal with 
the maximum at xo(a), and a;o(a) — >■ as a — > 0. 

The distributions ga{x) constitute basic ingredients of 
all theories of anomalous diffusion where they are em- 
ployed to produce solutions Pa{x,t) in space-time do- 
main of various forms of the Fokker-Planck equations 
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6,(5,1) 
&,(5,2) 
b,(5,3) 
6,(5,4) 



^^^(l) 

20ir S 
x/5-2''/^r(l) 

iov^r(^)i3 

3V5.3^/l»r(l) 

ior(A)r(T^)B 
^.3i/i0g-i/2^r(i) 



-1 

4^ 
207r A 

-v^.2-i/=r(§) 
iOv^r(^)A 

V5.37/ior(|) 

^or(A)r(il)^ 

6.27/105-3/2^r(f) 



V2r(j) 

167r 
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V^r(f) 

407r A 


-v^r(f) 

1207r S 


-V^.2l/5r(|) 

iO()^/^r(^)A 


v^-2-''/^r(l) 
i00v^r(^)s 


-3V5.3-Vi«r(#) 
25r(A)r(ii)A 


-7^.33/10r(i) 

75or(A)r(«)s 


14.23/105-5/2^r(f) 


ii.29/io5-V2^r(i) 



TABLE I. Coefficients of Eqs. © and ®; ^ = sin(7r/5) and B = sin(27r/5). 



along with their fractional generalizations (y, |2^- For 
instance in [22 Pa{x,t) is given as a convolution (called 
there inverse Levy transform) of ^ [—ga{t/s^/°')'\ with 
Pi{x,t) being a normalized solution of the ordinary 
Fokker-Planck equation, see Eq. ([Ij in [23l ■ The explicit 
forms of got (x) presented here will permit further devel- 
opment of this ambitious approach. 

The availability of giif. (x) makes it possible to fully de- 
scribe the long tail distributions of carrier transit times 
in amorphous materials such as As2Se3 and TNF-PVK. 
In fact in classic work [2J] the measured values of a for 
these two materials were a = 0.45 and a = 0.8 respec- 
tively, compare Fig. 6 of [2j|- These values were for a 
long time intractable theoretically. From now on, setting 
a = 9/20 and a = 4/5 in our Eqs. ([2]) and ^ directly 
provides the sought for framework for interpretation of 
these data. The appropriate distributions are presented 
as the curve II in Fig. ^, {a ^ 9/20) and the curve III 
in Fig. dSl), (a = 4/5). 



We believe that the exact forms of ga{x) obtained 
in this work, along with their asymptotics for a; — >■ oo 
and exact values of fractional moments, constitute a 
solid basis to extract a value of a best suited for an 
experimental situation at hand. Once it has been done, 
such description can be further 'fine-tuned' by choosing 
values of k and / which would optimise the choice of a. 
We hope that this approach will prove useful in practical 
applications. 

We thank Professor E. Barkai for kindly informing us 
about his results obtained in |22| . 
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FIG. 1. Comparison oi gi/k{x): the curves I, II and III corre- 
spond to l/k = 1/4 (s. Eq. (O), 1/3 and 1/2, respectively. 
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tively. Calculations were performed using Eqs. © and Q. 
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FIG. 3. Comparison of gi/f^ix): the curves I, II and III cor- 
respond to l/k = 2/5, 3/5 and 4/5, respectively. Calculations 
were performed using Eq. ([6]). 



to calculate Eq. ([2]): 

LevyDist := proc(k, 1, x) simplif y (convert ( 
sqrt(k * 1) * MeijerG([[] , [seq(jl/l, jl = 
0..1-1)]] , [[seq(j2/k, j2 = 0. .k-1)] . []] . 
I'l /(k'k * x'l)) /(x*(2*Pi)-((k-l)/2)), 
StandardFunctions)) : end: . 



